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Abstract. Using ISOPHOT maps at 100 and 200 yum and raster scans at 100, 120, 150 and 200 //m we have detected four 
unresolved far-infrared sources in the high latitude molecu l ar clo ud L 183. Two of the sources are identified with 1.3 mm 
continuum sources found by Ward-Thompson et al. (1994 , 2000) and are located near the temperature minimum and the 



coincident column density maximum of dust distribution. For these two sources, the ISO observations have enabled us to derive 
temperatures (~ 8.3 K) and masses (~ 1 .4 and 2.4 Mq). They are found to have masses greater than or comparable to their virial 
masses and are thus expected to undergo gravitational collapse. We classify them as pre-protostellar sources. The two new 
sources are good candidates for pre-protostellar sources or protostars within L 183. 

Key words. Stars: formation - ISM: clouds - dust, extinction ~ ISM: individual objects: Lynds 183 - ISM: individual objects: 
Lynds 134N - Infrared:ISM 



1. Introduction 

Low-mass stars are known to form within dark clouds. To study 
the initial conditions for their formation, it is important to probe 
the physical conditions of both molecular gas and cold dust 
in the deep interiors of such clouds. Dense cores have been 
identified inside dark clouds, some probably being in the stage 
of contraction to form a star. A pre-protostellar core is defined 
as the stage in which a gravitationally bound core has formed 
in a molecular cloud, but no central hydrostatic protostar exists 
yet within the core (see e.g. Ward-Thompson et al. 1994). Pre- 
protostellar cores are thought to be very cold. Thus, in many 
cases they escaped detection by IRAS, limited to A < 100 yum. 

The long wavelength and multi-filter capabilities of 
ISOPHOT aboard the Infrared Space Observatory (ISO) 
(Kessler et al. 1996| ), combined with its improved sensitivity 
and spatial resolution over IRAS, have been utilized in study- 
ing the far-IR emission of molecular clouds and pre-stellar and 
young embedded stellar objects within them. Analysis of their 
physical parameters is facilitated by using the spectral energy 
distributions (SEDs) obtained from ISOPHOT multi-filter pho- 
tometry. 



1.1. Lynds 183 

The dark cloud/large globule L 183, frequently cited also as 
L 134 North, is a prototypical dense cold molecular cloud. Its 
visual extinction is estimated to be ~ 17™ based on its 200 //m 



optical depth (Juvela et al. 2002). So far there is no evidence 



for associated newly born stars such as T Tauri stars or IRAS 
point sources. However, Martin & Kun (1996) have found a 



bona fide T Tauri star and Ha emission line star near L 183. 
These two stars are located outside our maps. Given its short 
distance of ~ 1 10 pc (Franco 1989 ) L 183 provides a good spa- 
tial resolution with the ISO FIR beam size (ISOPHOT's spatial 
resolution of 45" at 100 /im corresponds to 0.02 pc at the dis- 
tance of L 183). The location at high Galactic latitude (b = 36 
deg) minimizes contamination with unrelated ciiTus along the 
line of sight. The location off the Galactic plane at z = 65 pc 
implies that the impinging ultraviolet radiation field is strongly 
asymmetrical. 

A Digitized Sky Survey||]red plate ('Equatorial Red' survey 
(UK Schmidt) IIIaF + RG610) image of L 183 is shown in 
Fig. |l| In Fig. [l|b, the original image has been smoothed with 
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pressed digital form with the permission of these institutions 
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a FWHM=9"(nine image pixels) Gaussian kernel. The cloud 
has a dark core - bright rim structure, which is characteristic of 
clouds having high optical depth at visual wavelengths that are 
illuminated by the diffuse interstellar radiation field (Mattila 



Table 1. The properties of filters. A^^f is the reference wave- 
length, AA is the width, and dAiry is the diameter of the Airy 
disk 



1974; FitzGerald et al. 1976). The brightrim maximum occurs 
at a radius corresponding to extinction ss 1.5" - 2". 

Laureijs et al. (1991) observed the 60jum emission to de- 
cline relative to 100//m emission in a narrow transition layer 
around the cloud edge. They explained this behaviour by as- 
suming a separate grain component at 60 yum that undergoes a 
modification of properties in a transition layer. 



Laureijs et al. (1995) have used IRAS images, CO molecu- 
lar line observations and blue extinction values from star counts 
for a large-scale study of the L 134 cloud complex, includ- 
ing L 134, L 169 and L 183. Using '^CO observations, they 
found 18 clumps in the complex and derived their properties. 
The clumps follow clear size vs. linewidth and luminosity vs. 
size relationships. An analysis and discussion of the dust and 
molecular gas properties in L 183 based on the ISOPHOT 100 
and 200 //m emiss ion maps has recently been presented by 
Juvela et al. ( |2002| ). 

L 183 has been a favourite target for molecular line obser- 
vations which have revealed it as a rich source of molecular 
species (see e.g. Swade |l989a[ |l989b|; Dickens et al. |2000|). 



Lee et al. (2001) have classified L 183 as a strong infall candi- 
date. The optically thick CS(7=2-1) lines in the inner region of 
the cloud show a double peak with the blue component brighter 
than the red one, which is characteristic of inward motions. 
Similar profiles suggesting infall motions have been observed 
for HCO+(y=3-2) (Gregersen & Evans pOO0| ) and for CS(7=2- 
l)by Snellet al. ( |l982| ). 

Mapping of the millimetre and submillimetre continuum 
dust emission by Ward-Thompson et al. (1994, 1999) has re- 
vealed a small extended core in the center of L 183. This 
source has been detected at 800, 1100 and OOOjum, and has 
FWHM dimensions of 60" x40" (0.0 32 x 021 pc) at 800 jum. 
According to Ward-Thompson et al. ( 1994 ), the core is prob- 
ably pre-protostellar in nature, i.e. at an earlier stage than an 
accreting Class pr otostar (for the definition of a Class ob- 
ject see Andre et al. 1993| ). The submillimetre emission from 
the core is consistent with the dust being heated externally by 
the general interstellar radiation field. Recent observations by 
Ward-Thompson et al. (2000) have revealed that the contin- 
uum emission extends further to the south, and that there is 
another emission maximum located at a position about 1.5' 
south of the previously detected maximum, with a FWHM size 
of 120" X 60". Fr om su b-mm polarization observations Ward- 
Thompson et al. (2000) measured a magnetic field direction 
which is at an angle of 34° + 6° to the minor axis of the core, 
and found evidence for decreasing polarization at the highest 
continuum emission intensities. 



Ward-Thompson et al. (2002a) have mapped the core of 
L 183 at 90, 170 and 200 yum with ISOPHOT The core was 
not detected at 90 /im wavelength. The dust temperature de- 
rived from 170 and 200 //m data showed no temperature gradi- 
ent across the core. 

In this study, we search for pre- and protostellar ob- 
jects in the cloud. Our maps are much bigger than previous 
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(sub)millimetre continuum maps. The derived fluxes of the 
sources between 120-200 jum, located near the peak or in the 
Wien regime of the assumed blackbody radiation, are essen- 
tial for temperature and thus mass determination of the objects. 
Only if the masses of the objects are determined can we study 
their dynamical state. Furthermore, we study the relation of 
young (proto)stellar objects with the properties of the underly- 
ing dust in the cloud such as temperature and column density. 

In the context of L 183, the term 'pre-protostellar core' has 
been previously used to refer to the two known (sub)m m con- 
tinuu m sources within the cloud (Ward-Thompson et al. 1994 , 
2000), and also to the cloud core itself in which these sources 



are embedded (Ward-Thompson et al. 2002a ). Throughout this 
article, we use the term 'source' to refer to objects in L 183 
which are unresolved in our maps, such as the previously- 
known continuum sources. By the term 'core', we refer to the 
visually opaque cloud core which is larger than the sources 
within it and resolved by ISO observations. 

2. Observations and data reduction 

The observations were made with the ISOPHOT instrument 
aboard the ISO satellite, using the CI 00 and C200 detectors 
(Lemke et al. 1996). We have mapped ~ 30 x 30 arc min^ of 



the cloud at 100 and 200 /im. In addition, we have performed 
linear north-south raster scans at 100, 120, 150 and 200 //m 
through the cloud core with a two way path, the return path 
being shifted by 90" to the west (see Fig. ^for the raster scan 
positions and the orientation of the detector array). The raster 
scans were designed to cross the continuum source FIR 1 . All 
observations were made using the observing template PHT22 
in raster mode. Details of the detector/filter combinations and 
observational parameters are given in Tables [l| and |[ 

The data analysis was done using PIA|^ (ISOPHOT 
Interactive Analysis) V 8.1 (Gabriel et al. |1997| ). At the first 
processing level, the detector ramps were corrected for non- 
linearity of the detector response, glitches in ramps were re- 
moved by using the two-threshold glitch recognition method. 



- The ISOPHOT data presented in this paper were reduced us- 
ing PIA, which is a joint development by the ESA Astrophysics 
Division and the ISOPHOT Consortium (the ISOPHOT Consortium is 
led by the Max-Planck-Institute for Astronomy (MPIA), Heidelberg). 
Contributing ISOPHOT Consortium institutes are DIAS, RAL, AIP, 
MPIK, and MPIA 
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Fig. 1. Optical red Digitized Sky Survey (DSS) (a,b), 200 yum (c) and lOOjt/m (d) images of L 183. The four unresolved sources 
are marked as blue dots; for identification see Fig. ^ a) Digitized Sky Survey red plate. The image is scaled to best show the 
opaque core of the cloud, b) Digitized Sky Survey red plate. The image is scaled to best show details of faint surface brightness, 
c) False-colour image of 200 //m surface brightness. The colour scale is in units of MJy sr"\ d) False-colour image of 100 jum 
surface brightness. Zodiacal light of 4.0 MJy sr"' has been subtracted. The colour scale is in units of MJy sr"' 

Table 2. The parameters of the observations. TDT is the Target Dedicated Time number of the observation, P.A. is the position 
angle measured from North to East. All observations were made on January 19th 1997 
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and the ramps were fitted with 1st order polynomials. At subse- 
quent levels the signals were deglitched, reset interval correc- 
tion was applied, signals were linearized for the dependence 
of detector response on illumination, and orbital position- 
dependent dark currents were subtracted. 

For the flat-field correction, a statistical method was ap- 
pUed: the value of the reference pixel at each raster position 
was correlated with other pixels. Instead of comparing the ref- 
erence pixel only with pixels located at the same raster position, 
a mean of two pixels located symmetrically around the refer- 
ence pixel wass taken. In this way, it is possible to reduce the 
scatter caused by surface brightness gradients in the image. 

The primary intensity calibrators are the fine calibration 
source (FCS) measurements, which bracket the actual map or 
raster scan measurements. The ClOO detector suffers from re- 
sponse transients during the relatively short measurement time 
of the FCS. We have corrected for this effect by modelling the 
signal transients with the PIA signal drift interface, using the 
off'set exponential as the fitting function. The C200 detector re- 
acts faster, thus there is no such transient behaviour. 

We have compared our 100 /urn surface brightness cali- 
bration with COBE/DIRBE data. To match our data to the 
much coarser DIRBE resolution (0.7°x0.7°), we have used the 
IRAS/ISSA 100 jum map of the L 183 area as an intermediate 
step (for the method see Lehtinen et al. 2001). We have thus 
found a scaHng factor /io()^m(ISOPHOT)//ioo;,m(DIRBE) = 
1.25. Although this value is well within the combined uncer- 
tainty of the absolute accuracies cited for ISOPHOT (20%, see 
Klaas et al. pOOOb and for DIRBE (11%, see Hauser et al. 119981) 



it is larger than the scaling factor obtained for a sample of five 
mappings of different clouds where the ISOPHOT and DIRBE 
lOOyum surface brightnesses agree to within 15% (Lehtinen et 



al. 2001 ). We have thus rescaled our present ISOPHOT 100 yum 



map dividing it by 1.25 to correspond to the DIRBE calibra- 
tion. The North-South raster scan data at lOOyum are scaled to 
the map. 



A direct ISOPHOT vs. DIRBE comparison was done at 
200 fj.m by averaging the ISOPHOT surface brightness over our 
whole map area and comparing it with the nearest DIRBE pixel 
value. The resulting scaling factor was 0.90, i.e. well within the 
uncertainties of each instrument. No rescaling was applied to 
our ISOPHOT data at 200 jum, nor at 120 and 150/zm. 

3. Results 

3.1. Unresolved sources in L 183 

The ISOPHOT far-IR maps at 100/zm and 200 yum, together 
with optical images, are shown in Fig. |l]. We have searched 
for point-like objects by visual inspection. Neither one of the 
sub-mm continuum sources detected by Ward-Thompson et al. 
( |1994 |l999t |2000| ) is detected in the 100 yum map (see Fig. 0d). 
At 200 yum there is an emission maximum at the position of 
the sources (see Fig. They have an angular distance of ~1 
pixel, so they are not detected as separate sources in the 200 fim 
map. However, in the C200 North-South stripe scans the reso- 
lution is better due to overlapping pixels (see Fig. |^). A map 
made from the 200 yum stripe scan data clearly shows two emis- 
sion maxima whose positions correspond to positions of the 
continuum sources of Ward-Thompson et al. (see Fig. |^). In 
order to derive the total ffuxes of these sources, our sources 
FIR 1 and 2, we have modelled them by two-axial gaussian 
surfaces having positio n angles and FWHM sizes as given by 
Ward-Thompson et al. ( |1994 |2000| ). We have not made a 2-D 
map of the stripe scan data, but have treated them as a 1-D time- 
ordered data for each pixel. The background has been deter- 
mined by fitting the 1-D data with cubic splines, excluding the 
positions where the sources contribute to the ffux. The ffux seen 
by each pixel has been calculated by using our 2-D model and 
ISO beam profile. The ffuxes of the sources have been deter- 
mined by minimizing the difference between the modelled and 
observed ffuxes in a least squares sense. The model is shown in 
Fig. and the derived total ffuxes are given in Table |^. The 
accuracy of the derived ffuxes is mainly limited by our ability 
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15' 52" 30^ 00^ 5r30^ 00^ 
Right Ascension (1950.0) 

Fig. 3. Contour maps of 100 yum (a) and 200 yum (b) surface brightness. The detected unresolved sources are marked with asterisks. 
The additional raster scan positions (center of the detector array) observed at 100, 120, 150 and 200 ;um are indicated by dots. 
The orientation of the ClOO and C200 detector is shown for four raster scan positions. The sizes of the ClOO and C200 detector 
pixels are shown at the upper right corners, a) The contours are from 12 to 23 MJy sr ' in steps of 1 MJy sr ^ Zodiacal light of 
4.0 MJy sr"' has been subtracted, b) The contours are from 20 to 85 MJy sr ' in steps of 5 MJy sr ' . The chance location of FIR 3 
within the detector area is shown 

Table 3. Positions, sizes, and non colour-corrected total flux densities of unresolved sources in L 183. The 1-cr error values for 
the flux densities include an estimated 50% error for FIR 1 and 2 and 30% error for FIR 3 and 4 
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F,,(150yL(m) 
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[Jy] 
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15''51"'32':7 
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1 
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1 
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FIR 4 




15''51"48!4 
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<0.6 






4.7+1.6 





'Sizes and positions from Ward-Thompson et al. (1994, 2000) 



Table 4. The parameters for unresolved sources. Tj is the dust temperature. /Iq is the wavelength (corresponding to vq) where the 
optical depth is unity. The sub-mm luminosity Lsub-mm is the luminosity of the fitted blackbody function, integrated longward of 
350 fim. Lboi is the bolometric luminosity, N(H2) is the mean molecular hydrogen column density, and n(H2) is the mean number 
density assuming a constant density sphere with size equal to the geometric mean of the minor- and major-axis FWHM sizes. 
M is the total (gas plus dust) mass, and Mvh is the virial mass. The quoted errors are 1-cr errors based on Monte Carlo error 
estimation, including the estimated errors of fluxes only. Note that the parameters of FIR 3 have been calculated by assuming 
that its size is equal to the mean value of the sizes of FIR 1 and 2 
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Fig. 5. The spectral energy distributions of the unresolved sources. The 100, 120, 150 and 200 ;um flux density values are our 
ISOPHOT measurements, others are from the literature. All fluxes are total fluxes (see Sect. 3.2). The solid lines are modified 
blackbody fits (j6 = 2) to the data points, a) FIR 1. Derived temperature T=8.4K. The dotted lines are blackbodies which 
have 1 K higher and lower temperatures than the solid line, demonstrating the accuracy of the derived temperature, b) FIR 2. 
Derived temperature T=8.2K. c) FIR 3. Derived temperature T=12.8 K. d) FIR 4. The solid line is a modified blackbody curve 
(T=12.9K), with a temperature equal to the colour temperature derived from the 100 and 200 yum flux densities 



to subtract the strongly-varying background at the edge of the 
cloud. We therefore give an estimated error of ±50 % for the 
fluxes. 

A visual inspection of the 200 fim map further reveals two 
previously unknown unresolved sources in the cloud, sources 
FIR 3 and 4 (see Fig. [l|c). Neither of them is detected in the 
lOOjum map. FIR 3 is by chance located within the C200 de- 
tector array at one stripe scan position (see Fig. ||b), which 
enables us to also estimate its flux at 120 and 150yum. We take 
the flux of the pixel closest to the source, subtract from that the 
mean of the other three pixels, and correct the flux by the point 
spread function fraction for a point source centered on a pixel 
(Laureijs 1999 ). An error of 30% has been estimated for the 
fluxes. The 5.3 Jy 200 fim total flux derived this way is in good 
agreement with the 5.8 Jy flux derived from the map. Profiles 
through FIR 3 along image rows and columns in the 200 jum 
map are shown in Fig. ^ and d, respectively. 

Source FIR 4 is a prominent unresolved source on the 
200 /zm map (Figs. [l|c and On the 100 /zm map there is an 



extended brightness enhancement around the position of FIR 4 
(Fig. |l]d), where the 100 //m brightness reaches its maximum. 
In order to derive the lOOjum flux density, we have first con- 
volved the data to the 200 jum resolution. Since it is difficult 
to separate the source from the structured background, we con- 
sider the derived flux density as an upper limit. Profiles through 
FIR 4 along image rows and columns in the 200 /im map are 
shown in Fig. ^ and b, respectively. 

We have searched in the IRAS Point Source (PSC) and 
Faint Source Catalogs (FSC) for sources within our maps. We 
find two PSC sources, IRAS 15522-0258 and 15523-025 1, and 
one FSC source, F15507-0245. None of them are detected by 
us. The IRAS sources are located outside the cloud centre. The 
PSC sources only have 100 //m IRAS fluxes. In the ISO maps, 
there is no enhancement of lOOyum surface brightness at the 
position of the PSC sources. As the ISO and IRAS images at 
100//m appear morphologically similar, it seems that the two 
100 yum PSC sources are artefacts, and are a result of IRAS 
scanning dii^ection with respect to structures in the cloud. The 
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Fig. 2. (a) Map made from the 200 //m stripe scan data. The 
pixel size is 30". We note that this map is just to illustrate that 
the two cont inuum sources detected by Ward-Thompson et al. 
( 1994 2000 ) can be detected as separate sources in the stripe 
scan data. The details in this map are smaller than the actual 
resolution of the stripe scan data, (b) The model at 200 yum used 
to derive the fluxes of FIR 1 and FIR 2. The positions an d sizes 
of the sources are those given by Ward-Thompson et al. ( 1994 , 
200C). The inclined boxes show the locations of the 2x2 pixel 
C200 detector array in the North-South stripe scan. The con- 
tours of the continuum sources go from 0.8 to 6.4 in steps of 
0.8 Jy pixel ' . The FWHM of the ISO beam at 200 /urn is shown 
at lower left hand corner. The center position in both maps is 
the position of FIR 2 (see Table |]) 



FSC source has a 12 jjm detection only and is thus not expected 
to be detected at 100/zm. 

3.2. The spectral energy distributions 

At A ^ 100 /im, the emission arises from 'classical' large 
grains, which are at an equilibrium temperature within the radi- 
ation field. We have thus fitted the flux densities at these wave- 
lengths with a modified blackbody of the form 



= B,(rd)(l-exp(-T,))Q, 



(1) 



where By{Td) is the Planck function at the dust temperature T^, 
Ty is the optical depth which is assumed to vary with frequency 
as Ty = (v/vo'f, Qs is the solid angle of the emitting region, 
and vq is the frequency at which the optical depth is unity. We 
have assumed that the emissivity index /5 is 2. For the source 
size Qs of FIR 1 and 2 we have used the FWHM sizes given 



by Wai-d-Thompson et al. (1994; 200C). The size of FIR 3 is 



unknown, so we have assumed a size which is a mean value of 
the sizes of FIR 1 and 2. 



Fig. 4. Profiles along image rows and columns through the 
sources FIR 3 and 4 at 200 /im. The solid line is a fit to the 
background, the dotted line shows the fitted point source pro- 
file. The vertical line shows the fitted position of the source 



The 100 to 200 yum total fluxes measured with ISO are 
complemented by values at longer wavelengths from Ward- 
Thompson et al. (|1994|; |1999[) for FIR 1, and from Ward- 



Thompson et al. (2002b) for FIR 2. Ward-Thompson et al. 



(1994) give the total flux of FIR 1 only at 800 yum. From the 



given total flux and flux within a 18 "beam we derive their ra- 
tio of 10.4. Because the other fluxes at 450 and 1100//m are 
given for the same 18 "beam as the 800 yum flux, we can use 
the ratio 10.4 to derive total fluxes. This way we derive to- 
tal fluxes of <15.6 and 1.1 Jy at 45 and llOOyum, respec- 
tively. Ward-Thompson et al. ( 1999 ) have given for FIR 1 at 
1.3 mm a flux of 40mJy within 13"beam. Thus the total flux 
is (18/13)^ X 10.4 X 40mJy=0.8 Jy. The spectral energy distri- 
butions and fitted modified blackbody functions for FIR 1 , 2, 3 
and 4 are shown in Fig. |[ The resulting temperatures are listed 
together with other parameters in Table ^ The uncertainties 
of the fitted parameters have been estimated by Monte Carlo 
methods. 

For FIR 1, 2 and 3 the ISO flux densities, determined near 
the peak of the SEDs, accurately set the temperatures of the fit- 
ted blackbodies, as demonstrated in Fig. || The fitted tempera- 
tures are 8.4, 8.2 and 12.8 K, respectively. The derived colour 
temperature of FIR 4 has an upper limit of about 1 3 K. We 
note that the use of a single temperature is a simplification, 
and that the source sizes may be different at different wave- 



lengths. Ward-Thompson et al. ( 2002a i derived a temperature 
of 10±3 K within a 150"circular core area of the cloud, which 
is within errors equal to our value, 8.3 + 0.4 K. 
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3.3. The luminosities, column densities, number 
densities, and masses of unresolved sources 

The bolometric luminosities, Lboi, have been derived for 
all sources by integrating over the fitted blackbody curves 
(Table^. For the bolom etric lu minosity of the core of the cloud 
Ward-Thompson et al. ( 2002a ) derived 0.15 L0 by integrating 
the fluxes at 170 jum and 200 jum in a circular 150"aperture, 
combining these fluxes with mm/submm fluxes, and fitting the 
SED with a modified blackbody. The combined luminosity of 
IRS 1 and 2 derived by us, 0.09 Lq, is in excellent agreement 
with the core luminosity of Ward-Thompson et al. when tak- 
ing into account the different distances used by us (llOpc) and 
Ward-Thompson et al. (150 pc). 

The mean column densities through FIR 1 and 2 have been 
derived from the equation 



Nin2) = FJ(il,fimiiK,By(Ti)) 



(2) 



where Fy is the observed flux, fi - 2.33 is the mean molec- 
ular weight, mil is the mass of atomic hydrogen, and Ky the 
dust opacity per unit mass (gas-i-dust) column density (see 



e.g. Hildebrand |1983| ; Chini et al. |1987[ ). We have adopted 
fi300/jm = O.OOScnrg"' which is valid for pre-protostellar 
dense clumps (Motte et al. 1998 and references therein), as- 



suming a gas-to-dust ratio of 100. The value of Arnoo^m, and 
correspondingly the column de nsity, is uncertain by a fac- 
tor of 
N(H2) 



-2 (see e.g. Motte et al. 1998 ). We obtain the values 
^ 14 1022cni-2 and 8.1 lO^^^n-z for pjR 1 ^nd 2, re- 
spectively. These values are si milar to the values in cloud cores 
found by Motte et al. (1998 1 in p Ophiuchi cloud bas ed on 
millimeter continuum mapping, or by Bacmann et al. (200C) 
based on ISOCAM absorption survey. The difference by a fac- 
tor of about two in the column density of IRS 1 between us and 
Ward-Thompson et al. (1999) is due to the higher temperature 
(12.5 K) adopted by Ward-Thompson et al. 

Mean number densities n(H2), assuming a constant density, 
have been calculated using the formula «(H2) = (3/4)A^(H2)/^, 
where the cloud radius R for FIR 1 and 2 is assumed to be the 
geometric mean of the minor- and major-axis FWHM sizes. 

We are able to derive improved mass estimates because 
ISOPHOT flux values enable accurate determination of dust 
temperature. In order to estimate the total (gas plus dust) mass 
M of matter associated with dust we use the equation 



M : 



FyD^ 
KyBy(Td) 



(3) 



with D the distance. The fitted values of Aq, i.e. the wave- 
length where Tdust = 1 (see Table ^), suggest that for all 
sources the emission is optically thin at wavelengths ^100 fim. 
For FIR 1, we have used the 1300//m flux in Eq. (2), giving 
M X 1.4 M0. For sources FIR 2 and 3, we are using the 850 
and 200 fim fluxes, respectively. The value of k,i has been scaled 
from /fi3oo/im by using a A^^ dependence. The derived masses 
of FIR 2 and 3 are 2.4 and 0.2 Mq, respectively. 

The relative mass of the sources can be derived without the 
knowledge of or D. Only the dust temperature and flux at 
any optically thin wavelength are required. We have derived 
the relative mass between FIR 1 and 2 using the optically thin 



-02° 55' - 




15" 52"" 30^5^^ 00^ 5r45'' 30^^ 15-^ 
Right Ascension (1950.0) 



Fig. 6. Dust temperature as derived from the 100 and 200jt/m 
maps. The unresolved sources have been removed from the 
200 fj.m map before temperature calculation. The contours are 
from 11.4 to 14.2 in steps of 0.4 K. The detected sources are 
marked as stars. The temperature is not determinable outside 
the cloud due to the background subtraction. The size of the 
C200 detector is shown at the upper left corner 



1.3 mm fluxes from the blackbody fit. The derived mass ratio 
between FIR 1 and 2 is 0.6. 



3.4. Extended dust emission 

A detailed anal ysis of the extended dust emission is presented 
by Juvela et al. ( 2002 ). In this paper we will present the deriva- 
tion of the temperature and column density distribution of the 
dust. 

The deep red Schmidt image in Fig. |l] shows faint surface 
brightness details across the core, which are probably struc- 
tures in front of the dense core and are seen as bright regions 
due to scattered light. In the optical image, the eastern part of 
the cloud is covered by diffuse material. The dense part of the 
cloud which is delineated by the background stars extends un- 
till a « 15*'52™. Outside the diffuse material in the northern and 
western parts of the cloud the optical edge of the cloud closely 
follows the edge of 200 ;um surface brightness. The region of 
strong (^ 80MJysr"') 200 yum emission extends east behind 
the diffuse emission and follows the opaque cloud boundary. 
The diffuse emission is located in front of the cloud core. 

There is a major difference between the 100 and 200jt/m 
maps: the 200 /vm maximum around the sources FIR 1 and 
2 has no counterpart at 100 yum, i.e. there is a 200 /vm emis- 
sion excess relative to 100 /im. This could be due to a low dust 
temperature. Other explanations are that the dust emissivity at 
200 //m increases in the dense core due to grain coagulation 
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Fig. 7. Map of beam-averaged optical depth at 200 yum af- 
ter subtraction of the 4 point sources. The contours are from 
1x10"^ to 7x10"^ in steps of 1x10"^ 



(Cambresy et al. |2001t Stepnik et al. |2002| ), or that the two 
wavelengths trace different dust populations. 

In order to derive the colour temperature of dust, the fluxes 
at 100 and 200 /im have been fitted with a modified blackbody 
function of the form 



ly CC ^^By(Ti) 



(4) 



where /3 is the emissivity index and is the dust temperature. 
We convolve the lOOjum map to the 200 //m resolution, and 
subtract the background (determined at the upper right corner 
in Fig. 1^ and d). While fitting the modified blackbody, we it- 
eratively make colour corrections until the difference between 
subsequent temperatures is less than 0. 1 K. The temperature 
map (using y6 = 2) is shown in Fig.^ The minimum at the 
center of the cloud is ^ 11 . 1 K. 

In the optically thin regime, the optical depth is given by 
observed surface brightness 



(5) 



The map of beam-averaged optical depth at 200 //m is shown in 
Fig.0. The optical depth sharply decreases north of the sources 
FIR 1 and 2, while the decrease towards the south is smoother, 
similar to the morphology of the 200 jum map (Fig. ||c). The 
highest optical depth, T2oo/jm ~ 7 x 10"^, is ten times larger than 
the optical dept h thro ugh the center of the Thumbprint Nebula 
(Lehtinen et al. 1998), a globule without star formation. 



4. Discussion 

4.1. Are FIR 1 and 2 gravitationally bound objects ? 

By definition, a pre-protostellar object has to be gravitation- 
ally bound. Thus, the key question concerning the far-IR/sub- 



Table 5. The derived energies for FIR 1 and 2 in units of 10 J. 
The energy ratio is the ratio | Epot | /(Etherm + Eturb - Eext + 
O.SEmag) which has a value of 2 in the case of virial equilibrium. 
The magnetic energy E^ag is based on estimated upper limit for 
magnetic field strength B < lOO/zG 



Source 


Etherm 


Eturb 


Eext 




|Epot 1 


Energy 














ratio 


FIR 1 


15 


2.3 


0.35 


< 1.0 


91 


5.2 


FIR 2 


25 


3.4 


1.8 


< 5.4 


150 


5.1 



mm sources in L 183 is whether they are gravitationally bound, 
possibly already collapsing objects, or whether they are just 
some random density enhancements which will subsequently 
dissolve again. 

The improved mass estimates of FIR 1 and 2 enable us to 
make estimates of their dynamical state. A full treatment of 
virial equilibrium requires the consideration of gravity (fipot) 
and surface pressure (fiext) energies (compressing forces), and 
total kinetic (fikin) and magnetic (/imag) energies (supporting 
forces). The condition for virial equilibrium is then (McKee 
et al. 1993 and references therein) 



2(£'kin ^ext) "t" ^mag "t" ^pot 







(6) 



For the kinetic energy we have included the thermal energy 
of molecules and the turbulent energy. The thermal energy of 
molecules is 



-Etherm = 2^^^ 



(7) 



where N is the total number of molecules, and T is the kinetic 
temperature, and k is the Boltzmann constant. The turbulent 
energy is 

£.urb = ^Mal^ (8) 

where M is the cloud mass, and crtmb is the three-dimensional 
non-thermal (turbulent) velocity dispersion of the molecular 
gas, (Tturb - V^^^/(8 ln2) - kT/m, where AV is the observed 
linewidth (FWHM) and m is the mass of the molecule in con- 
sideration. We have estimated the velocity dispersion from NH3 
observations which have AV - 0.28 kms"' (Ungerechts 1980). 
For the kinetic temperature, we have used a value of lOK. We 
note that these values of velocity dispersion and kinetic temper- 
ature are values for the gas in which FIR 1 and 2 are embedded 
in. 

The external energy due to an external surface pressure act- 
ing on a core is approximately (Vallee 2000] ) 

- (i?^^)(?T^)'(T77^)' ^ 5 X 10^^ J (9) 
400cm ^ 0.6pc l.Okms ' 

where next and (Text are the density and velocity dispersion 
of the gas in which a core is embedded, and R is the core radius. 
We have used the values next = 3.0 x 10 cm -3 (Swade[l989§) 
and cText = 0.3 km s^' as giv en by i^CO(2-l) and C"^0(1- 
0) observations (Juvela et al. 2002). For R we have used the 



HWHM radii; 0.013 and 0.023 pc for FIR 1 and 2, respectively. 
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The total net magnetic energy is 
4nR^ 



3 Stt 



(10) 



which can be approximated by the formula (Schleuning 1998) 
B 



= (- 



R 



)^xlO^'^J 



(11) 



™^ MOOyuG' '^0.6pc' 
where B is the magnetic field strength. Recently, Crutcher 



& Troland ( 200C| ) have determined the line-of-sight magnetic 
field strength of ~ 1 1 jiG in the L 1544 core using OH Zeeman 
measurements with a 3' FWHM beam. Upper limits toward 
seven other dense cores gave Blos ^ 20 - 30yuG. Crutcher 
( 1999) has determined an upper limit Blos < 16 //G for L 183 
with a much larger (18') beam. Allowing for another fac- 
tor of four for the line-of-sight projection effects and for the 
higher density in the L 183 FIR 1 and FIR 2 cores we adopt 
B < 100 //G as an upper limit for the magnetic field strenght. 

The radial density pro file p{ r) of FIR 1 has been estimated 
by Ward-Thompson et al. ( 1994 ) based on the observed 800 
continuum flux density profile, and assuming that the core is 
isothermal; the density is p{r) oc 25 fj-Qjj^ j-jjg center to a 
distance of ~ 28" and p(r) cx furtherout. I n a sim ilar way, 
we have used Fig. 2 of Ward-Thompson et al. (2000) to derive 
the radial flux density profile for FIR 2. We find that a single 
power law for space density is not consistent with the observed 
radial flux density profile. Instead, two power laws with the 
same power law exponents as in the case of FIR 1 fit the data 
well. We have used these density profiles, with the addition of 
a constant-density central core having a radius of 1/10 th of the 
total radius. We have numerically calculated the gravitational 
potential energy Epot for FIR 1 and 2. For the sizes of FIR 1 and 
2, we have used the FWHM sizes given by Ward-Thompson 
et al. ( [1994| , |2000| ). 

The calculated energies are shown in Table |[ For both 
sources we have fitherm > ^turb > £^ext- The ratio | Epot I 
/(Etherm + Eturb - Egxt + O.SE^ag) h^s a value of two in virial 
equilibrium. For both FIR 1 and 2, the ratio has a value of about 
5. 

Neglecting the magnetic energy and external pressure, the 
virial mass can be written in the form 



= {(T^D)I{2G) 



(12) 



where D is the cloud diameter, G is the gravitational constant, 
and 0" is the three-dimensional velocity dispersion for a mean 
gas particle mass. 



IkT I Ay2 kT\ 
r \"h \81n(2) ~ 



(13) 



where m is the mean molecular mass, and m is the mass of 
the molecule used for observations. This form of virial equi- 
librium is exact for a polytrope of index n - 2 (Larson 1981 ). 
Using the values AV = 0.28 kms"' from NH3 observations and 
T=10K we obtain the virial masses 0.4 and 0.7 Mq for FIR 1 
and 2, respectively, which are about one third of the observed 
masses. Taking into account all uncertainties of the calculation, 
we conclude that the cores have masses which are higher than 
or comparable to their virial masses, and are thus gravitation- 
ally bound. 



4.2. Location of the sources relative to cold dust and 
gas 

The sources FIR 1 and 2 are located very close to the dust col- 
umn density maximum (see Fig. If we exclude the possibil- 
ity of a chance superposition, there is thus evidence that these 
sources were born at the very center of L 183. At the positions 
of FIR 3 and 4, the column density has decreased to about half 
of the maximum value. 

The NH3 emission is considered to be a good tracer of 
star-forming dense cores. The N H3(1,1 ) line in tensity maps of 
Ungerechts ( |198C| ) and Swade ( |19894 |l989b| ) show a north- 
south elongated structure with two maxima. The maxima are 
separated by about 3 arcminutes, and the sources FIR 1 and 
2 are located near the peak of the southern maximum. FIR 3 
is located at the very southern tip of the NH3 map of Swade 
(1989a). FIR 4 lies outside the existing NH3 and many other 
high-density tracer maps, but within the C'^O map of Juvela 
et al. ( |2002| ). 

In the C"*O(l-0) channel maps of Juvela et al. ( |2002| ) 
FIR 3 coincides with a clump which is most clearly sepa- 
rated at a velocity interval of 1.65-2.15kms"', while emis- 
sion around FIR 4 is strongest at velocities 2. 15-2.9 km s ^ 
Emission around FIR 1 and2is mostprominentbetween 2.15- 
2.9kms"'. At higher velocities, between 2. 9-3. 4km s"', the 
most intense emission is around FIR 1 and 2. It is thus pos- 
sible that FIR 3 is associated with a separate cloud component 
which is isolated both in spatial and velocity space. 



Juvela et al. (2002) present maps of H^^CO+(1-0) and 



DCO^(2-l) emission. FIR 3 and 4 are located outside the 
mapped areas. In particular the distribution of the cold gas 
tracer DCO^(2-l) looks very similar to the NH3 distribution; it 
has two maxima, the sources FIR 1 and 2 are located near the 
southern maximum, and the emission continues south towards 
FIR 3. The extent and structure of H^^CO^(l-O) emission are 
rather similar to DCO^(2-l); FIR 1 and 2 are located near its 
maximum and the emission continues towards FIR 3. 

4.3. Evidence for infall and outflow motions in L 183 



The detection of an optically thick molecular line (such as 
CS(2-1), HCO+(3-2)) with redshifted self-absorption, possi- 
bly with an absorption dip located at the same velocity as an 
optically thin line, is considered to be a strong evidence of in- 
ward motions. Recent observations of starless dense cores have 
shown that extended inward motions are common and may b e 
a necessary step in the star-formation process (Lee et al. 2001 ). 
The molecular line maps of L 183 by Lee et al. (2001) show 
CS(2-1) line profiles with infall asymmetry over an extended 
area of some 5'-6' around FIR 1 and 2, and N2H^(l-0 ) lines 
peaking at the absorption dip of CS(2-1) lines. Lee et al. (2001 ) 
have derived an infall radius of 3.8', which is about the same 
as the combined extension of FIR 1 and 2. 

One might wonder if the large-scale infall is connected 
with the build-up of the sources FIR 1 and 2, as they are the 



central sources of the core. Lee et al. (2001) have estimated 



the infall speeds and mass infall rates based on a simple two- 
layer radiative transfer model for clouds having infall asym- 
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Fig. 8. Molecular line profiles showing the infall asymmetry 
towards FIR 1 and 2. The CS(2-1) Une is from Lee et al. ( |1999| ), 
others are from us 



metries in their line profiles. They derive a mass infall rate for 
L 183 of ~ 1 - 2 X lO-^Moyr-i. With this rate, the build- 
up of the combined mass of FIR 1 and 2 would have taken 
~ 1.2 X 10^ - 2.3 X 10^ yr if the accretion rate is assumed to be 
constant. 

Fig . ^ shows the CS(2-1) spectrum from Lee et al. 
( |1999| ) (telescope FWHM«i27" at 86 GHz), and our DC O+(2- 
1), Hi3CO+(l-0) and €^^^8(2-1) spectra (Juvela et al. |2002| ) 
(FWHMsi57" at 86 GHz). Our spectra are at a single position 
nearest to the position of the spectrum from Lee et al. The lines 
are not directly comparable because the observations have dif- 
ferent angular resolutions. The CS(2-1) line is optically thick, 
while others are optically thin. The optically thin lines peak at 
the absorption dip of the CS(2-1) line at a velocity of about 
2.5 km s ', which confirms that the double-peaked profile of 
CS(2-1) line is not due to two cloud components at different 
velocities. In addition to the lines shown here, the optically 
thick HCO^(3-2) line shows a double-peaked profile with a 
brighter blue peak, and the optically thin H'^^CO^(3-2) line 



(Gregersen & Evans 2000) is peaking at the absorption dip 



at the same velocity as the other optically thin lines. This is 
consistent with the infall picture as derived from the CS(2-1) 
observations. 

The redshifted wing of the CS(2-1) line profile shown in 
Fig. ^may be an indication of outflow. Our '^CO(l-O) obser- 
vations (unpublished data) show many line profiles with either 
blue- or red-shifted wings or shoulders. However, the existence 
of blue- and red-shifted emission seems to be random, which 
could be due to two overlapping outflows, both from FIR 1 and 
2. 



Juvela et al. (2002) have studied the virial equilibrium of 
the L 183 cloud over a radius of ~ 10' (~ 0.3 pc) by consid- 
ering the kinetic and gravitational energies. They found that 
the cloud as a whole is approximately in virial equilibrium. If, 
however, the cloud is collapsing as suggested by asymmetric 
line profiles, the equilibrium has been unstable to gravitational 
collapse. This study of viiial equilibrium can not be apphed to 



the sources IRS 1 and 2 because their radii are much smaller, 
~ 0.02 pc. 

4.4. Nature of the sources 

Based on a submillimetre continuum stud y of 2 1 starless cores, 
including L 183, Ward-Thompson et al. ( 1994 ) concluded that 
the sources within cores are pre-protostellar rather than of 
accreting Class protostars. Firstly, the number of detected 
sources was much greater than expected from the estimated 
lifetime of Class objects. Secondly, the bolometric luminosi- 
ties of the sources were so low that if the luminosities were due 
solely to accretion, the ages of protostars would be too small 
compared to the number of sources detected. This conclusion 
was made on statistical grounds and some of the sources may 
indeed harbour a Class protostar. 

All the four sources in L 183 satisfy two defining properties 
of a Class protostar; 1) SED well fitted with a single modified 
blackbody (although the temperatures are lower than those of 
Class sources for which Td ~ 20 - 60 K (Andre et al. |200C| )), 
2) the ratio isub-mm/i'boi » 5 x lO"-'. However, a genuine pro- 
tostar requires further indirect evidence for the presence of a 
central protostellar object, such as a molecular outflow and/or a 
cm-wavelength continuum source. There is no confirmed out- 
flow in L 183 (see however Sect. 4.3). All confirmed Class 
objects, with one ex ception, show manifestation of an outflow 
(Andre et al. 2000 1. Harvey et al. (2002) have searched for 



3.6 cm continuum emission from possible protostars in L 183 
by using the NRAO Very Large Array. The half power diameter 
of the beam was 5.3' (i.e. effective field of view) but the spatial 
resolution of the VLA observations was ~ 11". No source was 
detected within the map of a size of about 8'square, centered 
near FIR 1 and 2 (the sources FIR 3 and 4 are outside or at the 
very edge of the map) with a 5-cr limiting flux level of about 
0.08 mJy. The seven Class candidate sources detected in Bok 
globules by Yun et al. (1996) have 3.6 cm continuum fluxes be- 
tween 0.2-5.3 mJy. It is thus most probable that FIR 1 and 2 
are pre-protostellar sources and not Class objects. 

The bolometric luminosities of FIR 1, 2 and 3 are very low, 
~0.05 Lq. In the list of confirmed Class protostars by Andre 
et al. ([2000|) the lowest Lbd is 0. 15 Lq. 



5. Summary and Conclusions 

The dark cloud L 183 has been mapped at far-infrared wave- 
lengths with the ISOPHOT instrument aboard ISO. Four un- 
resolved sources were detected in the cloud, called FIR 1, 2, 
3 and 4, of which FIR 1 and 2 were previously known. The 
SEDs of FIR 1 and 2 have been compiled by combining our 
ISO fluxes at 120, 150 and 200 ;um with longer wavelength 
fluxes from other studies. The SEDs are well fitted with a sin- 
gle modified blackbody with dust temperatures Tj ^ 8.3 K. 
The total masses (gas plus dust) of the sources are ~1.4M0 
and ~2.4Mq, and the bolometric luminosities are very low, 
-0.04 Lq and ~0.06Lq. Virial equilibrium consideration in- 
cluding magnetic, kinetic, and potential energy and external 
pressure shows that the sources have masses higher than or 
comparable to their virial masses. 
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FIR 3 is detected at 120, 150 and 200 jum. Its dust temper- 
ature is ~13K and bolometric luminosity is ~0.03Lq. FIR 4 
has been detected at 200 //m only. 

FIR 1 and 2 are probably gravitationally bound objects, 
rather pre-protostellar sources than Class protostars. The 
available data do not provide definitive conclusions on the na- 
ture of the objects. 

FIR 3 and 4 can be starless cold cores (either gravitation- 
ally bound or unbound), pre-protostellar cores (gravitationally 
bound) or protostars. Virial equilibrium considerations cannot 
be made due to a lack of accurate density structure and size 
information. Further (sub)mm and centimetre continuum ob- 
servations are required to study the nature of the sources. 
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